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We study single-hadron inclusive electroproduction in ep scattering at DESY 
HERA at next-to-leading order in the parton model of quantum chromodynam- 
ics endowed with non-perturbative fragmentation functions. Specifically, we con- 
sider charged-hadron production, with unspecified transverse momentum pt, in 
^ ' the Breit frame and D l± production as a function of px, and perform compar- 

l' isons with recent data from the HI Collaboration. 

i 1. Introduction 

In the framework of the parton model of quantum chromodynamics (QCD), 
. the inclusive production of single hadrons is described by means of frag- 

H ' mentation functions (FFs) D^(x,n). At leading order (LO), the value of 

D^(x, n) corresponds to the probability for the parton a produced at short 
distance to form a jet that includes the hadron h carrying the frac- 
tion x of the longitudinal momentum of a. Analogously, incoming hadrons 
and resolved photons are represented by (non-perturbative) parton density 
functions (PDFs) F a / h (x,fi). Unfortunately, it is not yet possible to cal- 
culate the FFs from first principles, in particular for hadrons with masses 
smaller than or comparable to the asymptotic scale parameter A. How- 
ever, given their x dependence at some energy scale the evolution with 
[i may be computed perturbatively in QCD using the time-like Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. Moreover, the factori- 
sation theorem guarantees that the D^(x,i.i) functions are independent of 
the process in which they have been determined and represent a universal 
property of h. This entitles us to transfer information on how a hadro- 
nises to h in a well-defined quantitative way from e + e~ annihilation, where 
the measurements are usually most precise, to other kinds of experiments, 
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such as photo-, lepto-, and hadroproduction. Recently, light-hadron FFs 
with complete quark flavour separation were determined 1 through a global 
fit to e + e~ data from LEP, PEP, and SLC thereby improving previous 
analyses. 2 ' 3 

In the following, we extend our previous report 4 on the electroproduc- 
tion, through deep-inelastic scattering (DIS), of n° mesons and charged 
hadrons with finite transverse momentum p^, m the j*p cm. frame at 
next-to-leading order (NLO) 5 by discussing charged hadrons with unspeci- 
fied values of p^,, including p* T = 0, and D*^ mesons with p* T > 0. 



2. Analytic Results 

At LO, inclusive hadron electroproduction proceeds through the Feynman 
diagram shown in Fig. 1(a), so that p* T — 0. At NLO, 6 virtual and real 
corrections, indicated in Figs. 1(b) and (c), respectively, contribute. In the 
latter case, p\. 1S integrated over. The NLO cross section is conveniently 
evaluated with the FORTRAN program CYCLOPS. 6 




The NLO analysis for the case that p^ > already at LO involves one 
more external parton leg and may be found in Refs. 5, 7. 
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3. Numerical Results 

3.1. Charged Hadrons in the Breit Frame 

HI 8 and ZEUS 9 measured the normalised Q distribution {1 / a^\s)da / dQ of 
charged hadrons in bins of x p = 2p Brelt /Q, where Q 2 = —q 2 is the virtuality 
of 7* and p Brolt is the projection of the three-momentum of h onto the flight 
direction of 7* in the Breit frame. In this frame, 7* is completely space- 
like, with four-momentum = (0,0, 0, — Q). This frame provides a clear 
separation of current and remnant jets and is especially appropriate for 
comparisons with inclusive hadron production by e+e - annihilation. 

In Fig. 2(a), preliminary HI data 8 are compared with NLO predictions 
evaluated with CTEQ6.1M 10 proton PDFs and AKK 1 FFs; the renormal- 
isation (r) and initial-state (i) and final-state (/) factorisation scales are 
taken to be /j, r = fa = fi r = £Q, where £ is varied between 1/2 and 2 about 
its default value 1 to estimate the unphysical-scale uncertainty. The PDF 
and FF uncertainties are assessed in Figs. 2(b) and (c) by switching to the 
MRST2004 11 PDFs and to the KKP 2 and K 3 FFs, respectively. 



3.2. D*± Mesons 

Among other things, HI 12 measured the p\, distribution da/dpip of D*^ 
mesons in the DIS range 2 < Q 2 < 100 GeV 2 and 0.05 < y < 0.7 with the 
acceptance cuts pr > 1.5 GeV and |?7| < 1.5 in the laboratory frame, where 
y is the relative lepton energy loss in the proton rest frame and r\ is the 
D*^ pseudorapidity. 

In Fig. 3, HI data 12 are compared with LO and NLO predictions eval- 
uated with CTEQ6 10 proton PDFs, BKK 13 FFs, and ^ 2 = = A- = 
£ Q 2 + (Pt) 2 f° r £ = 1- The theoretical uncertainty at NLO is esti- 
mated by varying £ between 1/2 and 2 about its default value 1. 



4. Conclusions 

We compared HI data on the electroproduction of charged hadrons 8 in 
the Breit frame and of D*^ 1 mesons 12 with p\, > with up-to-date NLO 
predictions. In the first case, we found reasonable agreement, except for the 
region of Q<30 GeV and x p >0.5, where the FFs are generally less well 
constrained by e + e~ data. In the second case, we found good agreement 
for p^> 1.25 GeV. This nicely supports the scaling violations in the FFs 
encoded via the DGLAP evolution as well as their universality predicted 
by the factorisation theorem. 
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Figure 2. The normalised Q distribution (l/<Tms)dcr /dQ of charged hadrons measured 
by HI 8 in bins of x p is compared with our NLO predictions estimating the theoretical 
uncertainties from the freedom of choice of (a) unphysical scales, (b) PDFs, and (c) FFs. 
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Figure 3. The pfp distribution da/dptp (in nb/GeV) of D*^ mesons measured by HI 12 
is compared with our LO and NLO predictions. 
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